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_. ■ Abstract 

P^ , The left tail of the implied volatility skew, coming from quotes on out-of-the-money put options, can 

'"^ ■ be thought to reflect the market's assessment of the risk of a huge drop in stock prices. We analyze how 

Cn , this market information can be integrated into the theoretical framework of convex monetary measures of 

risk. In particular, we make use of indifference pricing by dynamic convex risk measures, which are given 
as solutions of backward stochastic dififerential equations (BSDEs), to establish a link between these two 
^^ ■ approaches to risk measurement. We derive a characterization of the implied volatility in terms of the 

r^n , solution of a nonlinear PDE and provide a small time-to-maturity expansion and numerical solutions. 

This procedure allows to choose convex risk measures in a conveniently parametrized class, distorted 
^ 2 '. entropic dynamic risk measures, which we introduce here, such that the asymptotic volatility skew under 

^^ ■ indiff'erence pricing can be matched with the market skew. We demonstrate this in a cahbration exercise 

qh to market imphed volatility data. 
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1 Introduction 



X 



Risk measurement essentially conveys information about tails of distributions. However, that information 
is also contained in market prices of insurance securities that are contingent on a large (highly unlikely) 
5^ . downside, if we concede that those prices are mostly reflective of protection buyers' risk aversion. Examples 

are out-of-the-money put options that provide protection on large stock price drops, or senior tranches of 
CDOs that protect against the default risk of say 15 - 30% of investment grade US companies over a 5 year 
period 

A central regulatory and internal requirement in recent years, in the wake of a number of financial 
disasters and corporate scandals, has been that firms report a measure of the risk of their financial positions. 
The industry-standard risk measure, value-at-risk, is widely criticized for not being convex and thereby 
possibly penalizing diversification, and a number of natural problems arise: 
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1. How to construct risk measures with good properties. 

2. Computation of these risk measures under typical financial models. 

3. Choice: what is the "right" risk measure? 

The first issue has been extensively studied in the static case IIADEH991 IFS02II and recent developments in 
extending to dynamic risk measures with good time-consistency and/or recursive properties are discussed, 
e.g., in IIBK09I IKS05I IMZ091 IFSlll . However, concrete examples of dynamic, time-consistent convex 
risk measures are scarce, and they typically have to be defined abstractly, for example via the driver of a 
backward stochastic differential equation (BSDE) or as the limit of discrete time-consistent risk measures 
UStalOl . As a result, intuition is lost, and there is at present little understanding what the choice of driver 
says about the measure of risk. Or, to put it another way, how can the driver be constructed to be consistent 
with risk aversion reflected in the market? 

Let ^ be a bounded random variable representing a financial payoff whose value is known at some 
future time T < oo. A classical example of a convex risk measure, the entropic risk measure, is related to 
exponential utility: 

g(^) = -log{E[e-^^]), (1) 

where y > is a risk-aversion coefficient. When extending to dynamic risk measures gt(-) adapted to some 
filtration (;F/), a desirable property is (strong) time-consistency 

gA-gti^)) = Qsia 0<s<t<T. 

This flow property is important if Qt is used as a basis for a pricing system. The static entropic risk measure 
([T]) generalizes to 

gri^)^-log(E[e-y^\rt]). (2) 

The flow property follows simply from the tower property of conditional expectations. However, finding 
other directly-defined examples is not easy, and to have a reasonable class of choices, we need to resort to 
more abstract constructions. 

In a Brownian-based model, time-consistent dynamic risk measures can be built through BSDEs, as 
shown in IIBK09I IKSOSJ (compare also IIGia06ll ). extending the work of Peng IIPen04ll . That is, on a prob- 
ability space with a J-dimensional Brownian motion W that generates a filtration (T^t), the risk measure of 
the ?^r-measurable random variable ^ (e R for simplicity) is computed from the solution {Rt,Zt) € R x R'^ 
of the BSDE 

-dRt = g{t,Zt)dt-Z;dWt 
Rt = -^, 

where * denotes transpose. Here the driver g, which defines the risk measure, is Lipschitz and convex in z 
and satisfies g{t, 0) = 0. The solution is a process R e R that matches the terminal condition -^ on date T 
(when ^ is revealed and the risk is known), and a process Z e R^ that, roughly speaking, keeps the solution 
non-anticipating. Then Qt{^) := Rt defines a time-consistent dynamic convex risk measure. However, the 
possibility to offset risk by dynamically hedging in the market needs to be accounted for. Setting aside 
technicalities for the moment, this operation leads to a modification of the driver. 

The left tail of the implied volatility skew observed in equity markets is a reflection of the premium 
charged for out-of-the-money put options. The bulk of the skew reveals the heavy left tail in the risk-neutral 
density of the stock price 5r at expiration, but the very far left tail, where investor sentiment and crash-o- 
phobia takes over, could be interpreted as reveahng information about the representative market risk measure 



and its driver g, if we assume prices are consistent with this kind of pricing mechanism. The question then is 
to extract constraints on the driver from the observed tails of the skew, an inverse problem. In the application 
to equity options in Section [3l we assume the mid-market option prices reflect the premium a risk-averse 
buyer is willing to pay. We do not relate buyers' and sellers' prices to the bid-ask spread, since that is more 
likely related to the market maker's profit. 

To put our analysis into a broader framework, we observe that the underlying structural question is 
the inference of preference structures from observable data. The idea of using (at least in theory) observ- 
able consumption and investment streams to reveal the preference structure of a rational utility maximizing 
investor dates back to Samuelson in the 1940s and Black in the 1960s - for a recent overview on this "back- 
ward approach" to utility theory we refer to BCHQllj . The spirit of our presentation is a similar one, except 
we deal with dynamic risk measures rather than utility functions, and the observable data are not given as 
consumption and investment strategies but as readily available market implied volatilities. 

The main goal of the current article is to develop short-time asymptotics that can be used for the inverse 
problem of extracting information about the driver g from the observed skew. This could be used to construct 
an approximation to the driver and then to value more exotic derivatives in a way consistent with the risk 
measure. The information could also be used to quantify market perception of tail events, particularly when 
they depart from usual. Some studies have discussed the steepening of skew slopes in the run up to financial 
crises without a corresponding overall raise in volatility level. Infemng, fully or partially, a risk measure 
driver, could be used for detection of increased wariness of a crash. 

Berestycki et al. IIBBF04II presented short-time asymptotics for implied volatilities for no arbitrage 
pricing under a given risk-neutral measure in stochastic volatility models. Further work in this direction 
includes, among others, Wbbb 1 01 IFJ09 [ IFJLTOI and references therein. In Section |2l we extend this analysis 
to the nonlinear PDEs characterizing indifference pricing under dynamic convex risk measures. 

We find (Theorem l2. 121 ) that the zero-order term in the short-time approximation is the same as found in 
no-arbitrage pricing by IIBBF041 . The next order term is the solution of an inhomogeneous linear transport 
equation that sees only a particular slope of the partially Legendre-transformed driver, but is independent of 
the size of the options position (see equation (|27])). 

Section [3] illustrates the theoretical findings by focusing on a particular class of drivers, introducing 
distorted entropic convex dynamic risk measures. First we develop explicit calculations for the small time 
expansion in the Hull-White stochastic volatility model to illustrate the impact of the distortion parameter 
on the implied volatility skew. In Section 13.31 as a proof of concept, we perform a preliminary calibration 
exercise of the short time asymptotics to S&P 500 implied volatilities close to maturity. This allows to 
estimate the stochastic volatility model parameters from the liquid central part of the skew, and to recover 
the distortion parameter of our family of dynamic convex risk measures from far out-of-the-money put 
options. 

We illustrate the parameter impact for longer dated options in a numerical study (via the pricing PDF) of 
arctangent stochastic volatility driven by an Omstein-Uhlenbeck process. Section |4]contains the conclusions 
and Section [5] gives the more technical proofs omitted in the exposition. 

2 Heuristics and Statement of Results 

We consider a model of a financial market consisting of a risk- free bond bearing no interest and some stock 
following the stochastic volatility model on the filtered probability space {Q., T , {Tt), P) 



dSt =^i{Y,)Stdt + (T{Yt)StdW}, So = S; 

dY, - m{Y,) dt + a{Yt){p dW} + p' JW^), Fq = J, 



1 ^ „' A\\j1\ TZ _ .. (3) 



where W\ W^ are two independent Brownian motions generating (!Fr) and p' - yjl - p^. 



Assumption 2.1. VJe assume that: 

i) cr, a e C^^^ (R), where C^^^ (R) is the space of dijferentiable functions with locally Holder-continuous 
derivatives with Holder-exponent fi > 0; 

ii) both cr and a are bounded and bounded away from zero: 

0<o^<cr<o^<oo, and < a< a <a < oo; 

Hi) fj., m € C^^^^ (R), and \fi\ <J1 < oo. 

The pricing will be done via the indifference pricing mechanism for dynamic convex risk measures, 
which are introduced in the next subsection. 

2.1 Dynamic Convex Risk Measures, Indifference Pricing and BSDEs 

Definition 2.2. We call the family gt : U°{Q.,Tt,P) -^ U°{Q.,Tt,P), < t < T,a convex dynamic risk 
measure, if it satisfies for all t e [0, T] and all ^, ^', ^^ e L°°{Q., Tt, P) the following properties. 

/) Monotonicity: ^' > ^^ P-a..s implies Qti^^) < Qti^^)', 

ii) Cash invariance: gd^ + nit) = Qt(^) - nit for all m^ € L°°{Q., Tu P)', 

Hi) Convexity: gt{a^^ + (1 - a)^^) < agti^^) + (1 - a)gtif) for all a e [0, 1]; 

iv) Time-consistency gd^^) = gt{^^) implies gs{^^) = Qsi^^) for allO < s < t. 

We note that if the risk measure is additionally normalized, i.e. gt(0) = for all t e [0, T], then iv) 
is equivalent to the stronger property gsi-Qti^)) = Qsi^) for all < s < t IIKS05I Lemma 3.5]. The risk 
measure gd^) should be understood as the risk associated with the position ^ at time t. 

Ifgr is normalized, this is nothing else than the minimal capital requirement at time t to make the position 
riskless since gd^ + gd^)) = 0. In this static setting, the certainty equivalent price of a buyer of a derivative 
^ € L°°(D.,Tt,P) at time t is just the cash amount for which buying the derivative has equal risk to not 
buying it. 

In fact we are much more interested in the case where the buyer of the security is allowed to trade in 
the stock market to hedge her risk. In describing admissible strategies we follow the setting of continuous 
time arbitrage theory in the spirit of Delbaen-Schachermayer (for an overview, we refer to the monograph 
IIDS06II ). Denote therefore by 0^ the set of all admissible hedging strategies from time t onwards, i.e. all 
progressive processes such that 0, = and j 9s{pi{Ys) ds-\-o-{Ys) JW]) exists for all u e]t, T] and is uniformly 
bounded from below, and set 



-M 



%--\j ^.v(j"(n) ds + (r{Y,) dWl) : GeQ, 

The set of all superhedgeable payoffs is then given by Cr := (7C, - L°) n L°°, where L° denotes the set of all 
almost surely non-negative random variables. 

The residual risk at time t of the derivative ^ € L°°(n, Tj, P) after hedging is given by 

gr(^) := ess Mgd^ + h). (4) 

heC, 

Thus, assuming that the buyer's wealth at time t is x, her dynamic indifference price P,, which can be viewed 
as the certainty equivalent after optimal hedging in the underlying market, is given via gt(x + ^- Pi) - gt(x), 
whence, using cash invariance, 

P, = QtiO) - eti^l (5) 

4 



We note, while restricting ourselves to the buyer's indifference price, all our considerations ai^e easily adapt- 
able to the seller's indifference price by a simple change of signs of ^ and P, in dU. 

A convenient class of dynamic convex risk measures to which we will stick throughout this paper are 
defined from solutions of BSDEs. Assume that g : [0, T] X ^ x R^ ^ R is a f ® S{R^) predictable 
function which is continuous, convex, and quadratic (i.e. bounded in modulus by a quadratic function) in 
the R^-component. Next, let ^ e L°°(n, Tj, P) be a given bounded financial position. Then the BSDE 



R, = -€ + J g{t,aj,Zl,Z^)ds- J ZldWl-J 



Zl dW] (6) 



admits a unique y>-adapted solution {Rt, Z] , Z^), which defines a dynamic convex risk measure via ^r(X) := 
Rt IIBK091 Theorem 3.21]. 

The existence of a solution of the BSDE Q in this quadratic setting was first proved by Kobylanski 
HKobOOl Theorem 2.3], with some corrections to their arguments given by Nutz IINut07[ Theorem 3.6], 
whereas the uniqueness follows from the convexity of the driver as shown in IIBH08I Corollary 6]. From a 
financial perspective, the components Z\ Z^ of the "auxiliary" process Z can be interpreted as risk sources, 
describing the risk stemming from the traded asset and the volatility process respectively. 

2.2 Transformed BSDE under Hedging 

To assure the solvability of the BSDEs and PDEs that arise in our setting, we have to restrict the class of 
admissible drivers. Throughout, subscripts of functions indicate in the PDE context partial derivatives with 
the respect to the respective components. 

Definition 2.3. We call a function ^ : R^ — > R a strictly quadratic driver (normalized strictly quadratic 
driver) if it satisfies the following conditions i)-iii) (resp. o)-iii)): 

o) g(0,0)-0; 

/) g€C2-l(R2); 

") gzAZA fei , Z2) > for all (zi , za) e R^; 
///) there exist constants ci, C2 > such that 

^i(^ - (1 + ^2)) ^ «(zi'Z2) < C2(l + z? + z^) for all (zi,z2) e ^■ 
The normalization of the driver (condition o)) corresponds to the normalization of the risk measure. 

Remark 2.4. To ease the presentation, we only work with drivers that do not depend explicitly on time. 
While any dependence on Rt would destroy the cash invariance, it is not difficult to add an additional 
dependence of g on time (requiring some local Holder-continuity in the time component of g). The small 
time expansion Section 12.51 up to first order works in this case exactly as in the time independent case 
(evaluating the driver at t = T), and the higher order expansions can be adapted straightforwardly. For 
analysis of filtration-consistent, translation invariant nonlinear expectations, we refer to HCHYPOTj . 

In passing from the principal risk measure defined by g to the residual risk measure after hedging, as in 
(111), we will need the Fenchel-Legendre transform of g in its first component, namely 

g{^, Z2) := sup(^zi - gizi , zi)), ^ e R. (7) 

Lemma 2.5. Given that g is a (normalized) strictly quadratic driver, then the risk-adjusted driver g defined 
in (|7]) is also a (normalized) strictly quadratic driver 



Proof. To show g satisfies condition iii) of Definition I2.3[ we fix Z2 and treat tfie function as classical 
Fenchel-Legendre transform in one variable. Therefore it holds for proper, continuous convex functions 
/, g, that f < g imphes f > g and f = f, IIHULOlbl Proposition E. 1.3.1 and Corollary E. 1.3.6]. So the 
statement is proved by noting that 

sup(^zi - c(l + z? + zli) = c(^ - (1 + zl)) 

for any positive constant c. 

To show /) and //) in Definition 12.31 we note that condition ///) implies that g is 1-coercive in zi, i.e. 
g(zi,Z2)/\zi\ ^ oo as Zl — > ±oo for fixed zi- Now we can use the fact, that the Fenchel-Legendre transform 
of any 1-coercive, twice difFerentiable function with positive second derivative is itself 1-coercive and twice 
differentiable with positive second derivative, cf. BHULOlal Corollary X.4.2.10]. Thus it remains only to 
prove the differentiability of g with respect to Z2 which is a consequence of the differentiability properties of 
g: writing down the difference quotient and noting that the maximizer is differentiable, the positive second 
derivative with respect to the first component yields the existence of a finite limit. Finally g(0, 0) = follows 
from the definition if g(0, 0) = 0. n 

In other words, the class of strictly quadratic drivers is invariant under the convex conjugation in the first 
component and the class of normalized strictly quadratic drivers is an invariant subclass thereof. 

Remark 2.6. We note that it is important in our setting to stick to the theory of quadratic drivers, since if g 
would be a Lipschitz driver, g would be no more a proper function. This fact is easily seen, since from the 
Lipschitz condition it follows that 

g{Zl,Z2) < L(\ + yjzi+zl) < V2L(1 + Izil + |Z2|) 
for some constant L and hence 

g((,Z2) = sup(^zi - g(zi,Z2)) > sup(^zi - V2L(1 + Izil + IZ2I)) 

Zl Zl 

00 if\^\ > V2L 
-V2L(1 + |Z2|) if\^\< V2L. 

From now on we will assume that g is convex as a function on R^ and an admissible driver. Our next step 
is to describe the dynamic hedging risk in terms of BSDEs. These results are in essence due to Toussaint, 
IITou07[ Section 4.4.1]. Since his thesis is not easily available, we will nevertheless state the proofs here. It 
is convenient to introduce a notation for the Shaipe ratio: 

^Cv) := — — -■ (8) 

o-{y) 

Proposition 2.7. The risk of the financial position ^ e L"'{D,,T^r,P) under hedging is Qt{^) = Rt where 
Rf is given via the unique solution of the BSDE 



Rf = -^-J ZlA(Y,)+g{-A{Y,,),Z^)ds- J ZldWl-J Z^dWl 



(9) 



Moreover, gt is itself a dynamic convex risk measure. 



Proof. It follows from the work of Kloppel and Schweizer IIKS051 Theorem 7.17] that the risk is given via 
the BSDE 

Rf = -^ + ^^ ~g(Z],Zl)ds- ^^ Z]dWl- ^^ Z^dWl 

where g is given by the infimal convolution 

gizuzi) ■- inf(^(zi + cT{Y,)v,Z2)+KYt)v). (10) 

To be precise, besides the differences in sign between our convex risk measures and their monetary concave 
utility functionals, our L°° framework is in line with the main part of their paper where they work in L°° . 
However, the result IIKS05I Theorem 7. 17] is stated in the framework of L^-BSDEs with Lipschitz drivers. 
Their detour to L^ was due to their consideration that this is the natural framework for BSDEs. We have 
motivated that we have to work with quadratic drivers, for which there is yet no L^-theory, but it is straight- 
forward (though tedious) to check that their result (flOl ) adapts to our setting due to the regularity enforced 
by the admissibility conditions in Definition l2.3l 
Next, we rewrite the infimum in (ITOl) to get 



g(zi,Z2) = ini(giu,Z2) - (zi - u)A(Y,)) = -ziA(Y,) - g{-A(Y,),Z2). 

Finally, the uniqueness of the solution of the BSDE ^ follows again from IIBHOSi Corollary 6] using the 
convexity of the driver, which is implied by the fact that g is concave in the second component. Moreover, 
this entails also that ^ is a dynamic convex risk measure. n 

Remark 2.8. 1. We remark that in view of equation (O of the above proposition, the notion of admis- 
sibility could be slightly extended: it is possible to replace the lower bound in Definition I2.il Hi) 
by 



Cl(/(Zl)-(l+Z^))<g(Zl,Z2) 



for an arbitrary real-valued, convex and l-coercive function f. This is enough to get existence and 
uniqueness of equation ^, however it would clearly destroy the nice invariance property of Lemma 
\2.5\ and we do not adopt it in the following. 

2. As in ^ and the rest of the paper, the first argument of g is the negative of the Sharpe ratio. If we 
knew a priori that the Sharpe ratio was small relative to the minimal Lipschitz constant L in Remark 
\2.6\ then we could allow BSDEs with Lipschitz drivers. However we choose to put restrictions on the 
driver, such as in Definition \2.3\ iii ). rather than on the range of the model parameters. 

2.3 Indifference Valuation of European Claims 

From the formula (|5]l for the indifference price P, and Proposition 12.71 we have that 

P, = Rf'-Rf. (11) 

From now on we will restrict ourselves to particular bounded payoffs, namely European put options with 
strike price K and maturity date T: ^ = {K - StY . Moreover, for the further treatment the substitutions 

x:=\og{SIK), T-T-t, (12) 

will be convenient and we introduce the following notation. Denote by Lj the layer [0, T] x M? and by 
2ro,r the open cylinder above the disk B{m, r) with midpoint m, radius r and height < tq < T: QTo,r '■= 



]0, To[xS(m, r). Since the location of the midpoint (once fixed) will play no further role, we skip it in the 
notation. 

The following theorem characterizes the indifference price of a European put with respect to the dynamic 
convex risk measure with driver g under the stochastic volatility model ^. 

Theorem 2.9. The buyer's indifference price of the European put option is given as 

P(T,x,y) = u{T,y)-u(T,x,y) (13) 

where u € C'^^'^'^"^^((2r,r) ^ C{LT)for every cylinder Qj-,- is the solution of the semi-linear parabolic PDE 

{ -Ur + Lu= ■^g{-A(y),p'Ka(y)Uy); 

\ u{0,x,y)^ -(1-rT, ^ ' 

with operator L given by 

1 , 

Lu - Mu - -cr {y)Ux + {m(y) - pa(y)A{y))uy 

1 2 1 2 

Mu = -cr (y)urx + po-{y)a(y)Uxy + -a (y)Uyy (15) 

and u denotes the (x-independent) solution of (1141 ). with altered initial condition u{0, y) = 0. Moreover there 
exists a solution of the Cauchy problem (1141 ) (as well as that with the altered initial condition) which is the 
unique classical solution that is bounded in Lj together with its derivatives. 

Proof. By Ladyzhenskaya et al. IILSU67[ Theorem V.8. 1 and Remark V.8. 1], there exists a classical solution 
V € C^+'^/2,2+/?(-g^ ) ^ (-(-^^) jQ jj^g semilinear parabolic PDE 

1 1 

Vt + ■zcr^'iy)S'^vss + -zaiyfvyy + po-{y)a{y)Svsy + {m{y) - pa(y)A(y))vy = g{-A(y),p'a(y)Vy), (16) 

with terminal condition v{T,S,y) = -{K - 5)^ which is the unique classical solution to this PDE that is 
bounded in Lj- together with its (first and second order) derivatives. Applying Ito's formula to v{t,S t, Yt) 
and defining 



7' 

Rt 



= (T{Yt)StvAt,St, Yt) + pa(Y,)vy{t,St, Y^) 
= p'a{Yt)Vy{t,Sr,Yt), 
= v{t,St,Yr) 



shows that {RuZ},Zj) solves the BSDE ^ for {Rf,Z},Z^) with ^ = {K-STT,and therefore we identify 
R\ - v{t,St, Yf). The transformation (fT2l ). together with u{t, x,y) - v{t, s,y)/K leads to the Cauchy problem 
(fT4l ) for u. Finally, taking zero terminal condition for the PDE ( fT6l ). and calling the solution v(t,y) leads to 
Rf = v{t, Yf). Therefore the indifference price in (fTTI ) is given by Pf = v{t, Yt) - v{t,St, Yt), which, in 
transformed notation, leads to (fT3l) . n 



Here u is the value function of the holder of the put option, and u is related to the investment (or Merton) 
problem with trading only in the underlying stock and money market account. The nonlinearity in the PDE 
(fT4l ) is in its "Greek" Uy, that is, the Vega, and enters through the Legendre transform of the driver g in 
its first variable. For the familiar entropic risk measure, g(zi,Z2) = yiz^ + Zj)/'^, where y > is a risk- 
aversion parameter, we have g{^, Z2) = i^^/y - yZj)/'^- In this case, the nonlinearity is as m^ (see for example 
llBK05l[SZ05l ). 



2.4 Implied Volatility PDE 

Our main goal is to establish an asymptotic expansion of the indifference price implied volatility in the limit 
of short time-to-maturity. To do so, we now adapt the approach of Berestycki et al. IIBBF04II to establish 
a PDE satisfied by the Black-Scholes volatility /(t, x,y) implied by the indiff"erence pricing. Therefore we 
note that in the Black-Scholes model with unit volatility, the no arbitrage pricing PDE is given by 

UiO,x) =(l-rT, 
which can be represented explicitly as 

where <1> is the cumulative density function of the standard normal distribution. Using the scaling properties 
of the Black-Scholes put price, the indifference pricing implied volatility I{T,x,y) is hence given by the 
equation 

P{T,x,y) = Uil\T,x,y)T,x). (17) 

Before we derive a PDE for the implied volatility, we give some a priori bounds on indifference prices and 
their associated implied volatilities. 

Proposition 2.10. Denote by P^^{t, x; cr) the Black-Scholes price of the put calculated with constant volatil- 
ity a. Then 

P^^ (t, x; a) < P{t, X, y) < P^^ (t, x;W) (1 8) 

and 

g^<I{T,x,y)<o', (19) 

where a_ and a are the volatility bounds in Assumption \2.1\ 

The proof is given in Section |5] 

To derive the PDE for the implied volatility, we plug (ITT] ) into the equation (IT3] ) and get after some 
calculations - detailed in Section [5] - that the implied volatility / is subject to the nonlinear- degenerate 
parabolic PDE 

. 72x jKAj r2^-^ 1 2r2^r r^n ^ .. g{-^(y)'P'Ka{y)iiy)-g{-A(y),p'Ka(y)uy) 
- (rnr + tIMI + I'g- - -T^rGI + Tq{y)Uy = 2tII,. -— -, (20) 

7 4 ■ K{Uy - Uy) 

with the square field type operator 

gi ■- Mf - 2IMI (21) 

and 

q(y) := po-(y)aiy) -h 2m(y) - 2pa(y)A(y). (22) 

To motivate an initial condition for the Cauchy problem, we send formally r to zero in (l20l ). As the 
"Vega" V = Uy- Uy tends also to zero as r i (this is shown in Lemma [5^ . we observe that the quotient on 
the right side of (I20l) 

g{-A(y),p'Ka(y)Uy)-g{-A(y),p'Ka(y)Uy) -,, ^ n^ 
-— ■■ > p a(y)g,^{-A{y),0), 

K(Uy - Uy) 

which is bounded by the definition of strictly quadratic drivers and Lemma 1231 Dividing by 7^, this leads to 
the formal limit equation 

g 1. (23) 

mx,y) 



Remark 2.11. Our Cauchy problem is similar to that derived for no arbitrage pricing implied volatilities 
in l[BBF04\l . where they have the same equation (1201 ). but i) without the last term on the left side (which 
here is due to the change in measure from physical to a risk-neutral one); and ii) without the right side term 
(which here is due to the dynamic convex risk measure used for indifference pricing). However, our initial 
condition (I23I ). which does not depend on g and the drift of the stochastic volatility model, is exactly the 
same as theirs. 

Now we turn this heuristic argument into a precise statement. 

Theorem 2.12. The implied volatility function I{t, x,y) generated by the indifference pricing mechanism is 
the unique solution I e C^^^'^'^^^iQj^r) ^ C{Lt) to the following nonlinear parabolic Cauchy problem 

/ 72^ ,.,, ,2^-^ 1 2,2^, , xrr g{-A{y),p'Ka{y)Uy)-g{-A{y),p'Ka{y)u;) 
- (THr + tIMI + rg- - -rT-'rgi + Tq(y)IIy - iTlIy -— (24) 

7 4 K{Uy - Uy) 

where At, Q, and q are given by ( 1151 ). (1211 ). and (1221) respectively. The initial condition is given as 7(0, x, y) - 
x/il/ix,y) where >p is the unique viscosity solution of the eikonal equation 

&^ =1; 

<A(0,y) =0; (25) 

\p{x, y) > for X > 0. 

The proof is given in Section |5l 

It is worthwhile to note that indifference prices are not linear. Indeed, buying double the amount of 
securities will not lead to twice the price. In this way also the volatility implied by indifference prices is 
quantity-dependent as a consequence of the appearance of Uy and Uy on the right side of equation (l24l ). 
However, the nonlinearity in quantity is not observed in the zeroth- and first-order small time-to-maturity 
approximation as we show in the following subsection. Moreover, deriving the PDE for the indifference 
price implied volatility for buying n put options results in the same Cauchy problem as in Theorem 12. 12[ 
where one has only to replace K by nK in every appearance in equation (l24l ). 

2.5 Small-Time Expansion 

In the short time limit the implied volatility under indifference pricing is equal to the usual one as calculated 
by Berestycki e.a. IIBBF04II . as the initial conditions are the same (see Remark l2.11l) . The subtleties of the 
indifference pricing appear only away from maturity. This can be seen by methods of a formal small-time 
expansions. We make the Ansatz of an asymptotic expansion of the implied volatility: 

7(r, x,y) = I°{x,y){\ + Tl\x,y) + 0{t^)). (26) 



As seen above, the term I^\x,y) = 1(0, x,y) is given via solution of the eikonal equation (1251 ). To find the 
PDE for 7', we plug the expansion (l26l) into the equation (|24l) and compare the first order terms for t — > 0. 
This leads to the inhomogeneous linear transport equation 

27'(x,3;) + '^((r\y)i^,+pcr(y)a(y)i/,y)ll{x,y) + '^(pcr(y)a(y)ifr, + aHy)i^y)ll,(x,y) - ^M^ 

= - {q(y) - p'a(y)g,,i-A(y), 0))^ (27) 



where iff is again the solution of the eikonal equation (1251) . 
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It is important to observe that the dependence of this first order approximation on the risk measure (via 
its driver g) is given merely by the evaluation at Z2 - of the derivative of its Fenchel-Legendre transform 
g with respect to the second component. 

Comparing our PDE to the analogous equation in the arbitrage-pricing setting of IIBBF04II (who, how- 
ever, prescribe only the methodology in general and make explicit calculations just in one example), we note 
the presence of two additional terms in (|27] ) - one due to the change to a risk-neutral probability measure 
and the second due to indifference pricing with a dynamic convex risk measure. 

Furthermore in many we can obtain an interior boundary condition for the PDE at jc = by sending 
X formally to zero in (|27]) . since by i/^ ^ and ipx it follows that the first order coefficients vanish. We 
will show this specific procedure on a concrete example in Section 13.21 Imposing higher regularity on the 
coefficients and on g one can obtain also higher order terms in the expansion of the implied volatility. This 
is done by using Taylor expansions of g (in the second component) and Uy and Uy (in r). 

3 Examples & Computations 

In this Section, we introduce a family of dynamic risk measures within which to present the effect of risk 
aversion on implied volatilities, first using the asymptotic approximation in the Hull-White stochastic volatil- 
ity model, and later using a numerical solution of the quasilinear option pricing PDE. 

3.1 Distorted Entropic Dynamic Convex Risk Measures 

To study the impact of the driver on the implied volatility, we will now turn to a nicely parametrized family 
of risk measures. Therefore we define the following class of drivers, generating distorted forms of the 
entropic risk measure: 

2 

g^'-^zuzi) - \{z\ + zl) + mz\Z2 + "^zl = ^((zi + rizif + zl). (28) 

It is clear, that in the case 77 = ffiis is ffie driver connected to the classical entropic risk measure, whereas t] 
can be seen as a parameter which describes in which way volatility risk increases also the risk coming from 
the tradable asset. As we will see later in Section [l!4l 77 effectively plays the role of a volatility risk premium. 
In the case of the usual entropic risk measure the driver describes a circle whose radius is governed by the 
parameter y. In the distorted case it is now an ellipse where 77 determines additionally the eccentricity. 
Turning to ffie Fenchel-Legendre transform, we have 

1.2 _ r. 

2y^ r 
Plugging this into (l24l ). we see that the right hand side now reads 

Tllyilrip A.{y)a{y) - yKp'^a(yf{Uy + Uyfj 

and we remark in particular that y scales with K and hence also with the number of securities bought (as 
mentioned at the end of Section 12.4b . In particular we see again that the term appearing in the first order 
approximation of the implied volatility, g1^i-^{y), 0) = T]A(y), is independent of y. 

3.2 Short-Time Asymptotics for the Hull- White Model 

In the following we look at an example which is an adaption of Example 6.1/6.3 of IIBBF041 . Let ffie 
stochastic volatility model be given as the Hull- White model 

dSt ^KYt)Stdt + Y,S,dW}, So = s; 

dY, =KY,dWf, Yo = y. ^^ 
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f'H^,Z2) = —^^-^zl-rj^Z2. (29) 



for two independent Brownian motions W', W^. Obviously the model does not fall in the class considered 
above because the volatility cr(Y) = 7 is a geometric Brownian motion that is not bounded above or away 
from zero. Nevertheless, by a change of variables we will see that the results hold. 

Writing down the pricing PDE in the case of the distorted entropic risk measure (1281 ). we get 

1/2 2 2 ^ 1 2 1 y"Cv)^ yKK^ 22 , . 

-Ur+^yy "v.v + K y Uyy) - -y M.v ^ Y^ 2 ~y ">' "*" mi^W^y, 

and one sees that by the time change r i-> Ty^ (the boundary j = is not hit when we start with jo > given 
that Yf is a geometric Brownian motion) and setting JKy) :- fi(y)/y^ the equation becomes 

1 2x1 1 2 yKK^ 2 

-Mr + ^y'^xx + K Uyy) - -U^ ^ Y^^'^^ 2~">' "*" ^'^f^'^^'^y 

This equation has a solution (again by IILSU67I Theorem V.8.1 and Remark V.8.1]), at least in the case that 
fi is locally y6-Holder continuous (which in turn implies that ij.(y) - 0{y^) as y ^ 0) . 

In absence of global bounds on the volatility we are not able to derive results about existence and 
uniqueness of solutions for the PDE (l24l) for the implied volatility. Nevertheless we can postulate the small- 
time expansion to get 

' «A? + ^¥v =1// 

<A(0,y) =0 (31) 

\p{x, y) > for X > 

as the PDE characterizing it's zeroth oder term and 




12 ' '' 2.V 

for the first. The interior boundary condition for /' at x = follows from the formal asymptotics as x — > 0. 
As derived in IIBBF04I . the zeroth order term of the expansion is given via the solution of ( [3T] ). 



^(x,y)^lln(^+ Jl + ^). (33) 

as /°(x,y) - x/ilf(x,y), whereas for (l32l ) we can guarantee only a solution in the case where /i(y) - 0(y^) 
since I^/I^ = -ifry/ifr = o{l/y) as y ^ 0. Obviously this means practically that we need an extreme drift in 
the Hull-White model to compensate the very volatile volatility process. However, setting e.g. fi(y) = jiy^ 
for some constant //, we are able to solve the PDE ( [32l ) explicitly by the method of characteristics to get 

In Figure [H we rely on the parameter set 

^ = 6; K = l; y = 03; t = 0.1. 

Whereas the parameter y does not appeal" in the first order approximation, the distortion parameter 77 has 

a double effect. On the one hand side it shifts the smile at the money a small amount, on the other hand it 

2 9 

changes more significantly the wing behavior of the smile, adding to the asymptotics the term 77^ . "j , 2 
(since x^/i/r^ ~ ^^x^/(ln \x\f- as x ^ ±00). This changes the whole wing behavior, since 7° ~ ^1^ ^rid 
/' — > for 77 = as X ^ ±00. Of course 77 = corresponds to the first order term of martingale pricing as 
IIBBF041 . Positive 77 (hence a positive impact of the volatility risk on the risk of the traded asset) increases 
the implied volatility and steepens the wings. 
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Figure 1: Implied volatility in terms of log-moneyness log(.K'/5'o) for the Hull-White model: zeroth and first 
order approximation in dependence of 77. 



3.3 Calibration Exercise 

We perform a simple calibration exercise using S&P 500 implied volatilities to show how the approximation 
(l26l ) could be used to recover some information about the market's risk measure. We work with the short 
time approximation in the uncorrelated Hull-White model (l30l) of Section 13.21 and within the family of 
distorted dynamic convex risk measure (l28l) . 

We suppose that the near-the-money option prices ai^e given by expectations under an equivalent martin- 
gale measure, specifically the minimal martingale measure. This corresponds to setting the distortion 77 = 0, 
and so the short-time asymptotic approximation for implied volatilities / is 

where \li{x,y) is given by (1331 ) and l'^{x,y) by (l34l) with 77 = 0. This gives 



I^Im-- 



Hx,y) 



1 + 



ilj\x,y) 



In 



y / K X 



.2 ^2 \ 4 



where we see that the approximation Im depends only on the parameter k and the current volatility level y. 
Then having fit the liquid implied volatilities to recover estimates of k and y, we can write the short time 
approximation with distortion as 



I X Im + iinf)- 



TX 



(35) 



2^\x,yy 

from which we can estimate the combination iirj. 

In Figure [21 we show the result of a fit to S&P 500 implied volatilities from options with nine days to 
maturity on October 7th, 2010. We comment that the Hull-White model is not the best choice of stochas- 
tic volatility model, and our intention is primarily to demonstrate the procedure. In particular, because 
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our asymptotic formulas are for the uncorrelated model, which always has a smile with minimum implied 
volatility at the money, we restrict the fit to options with negative log-moneyness, which is mainly out-of- 
the-money put options. The first part of the least-squares fitting to options with log-moneyness between 
-0.06 and gives k ~ 6.6, and 3^ « 0.18. Then fitting to (1351 ) the options with log-moneyness less than 
-0.06, which we view as less liquid and more reflective of the market's risk measure, gives the estimate 
HT] w 35. We note that the n in this model is expected to be rather large because the composite drift of the 
stock is //y^, and y « 1. So the estimate we get from the data of ^^77 is not unreasonable. In particular it 
reveals a positive distortion coefficient 77. 




■Risk measure fit to far oul-of-the-money puts 



■Fit to iiquid option impiied vois 




-0.06 -0.05 -0.04 -0.03 -0.02 -0.01 

Log-Moneyness 



Figure 2: Fit of short-time asymptotics to S&P 500 implied volatility data 9 days from maturity. 



3.4 Numerical Study 

We consider the buyer's indifference price of one European put option with respect to the family of distorted 
entropic risk measures defined by (|28] ). We work within the stochastic volatility model Q and, for the 
numerical solution, we return to the primitive variables (f, S,y). Denote by X.y the generator of the Markov 
process F: 



■Cy = 2«0') TT + "^Cv)— ' 



d_ 
'dy 



and by Xs,,, the generator of (5, F): 



1 



Ls,y = £y + ^^iyfS^^^ +pa{y)(r{y)S 



.id' 
55^ 



52 d 

dSd-y^^^W 
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From (fTTl) . the buyer's indifference price of a put option with strike K and expiration date T at time t < T 
when 5, = S and Yt - y, is given by: 

P{t,S,y) = ip{t,S,y)-^o{t,y), 

where ip solves 

ft + \Ls,y - pa{y)^(y)-T- V ^ -f (-^(y), -p'a(y)(py), (36) 

ip{T,S,y) - {K-S)\ 

and </'o(f^3') solves the same PDE without the 5 -derivatives and with zero terminal condition. Note that 
(fi = -V, where v was the solution to the PDE problem in (IT6l) . and (fo = -v which was introduced in the 
proof of Theorem l2.9l 

As g is given by ( [291 ). we can re-write (l36l ) as 

/ , d\ A^(y) 1 , , , 

ft + Us,, - (P + vpMyUiy)-\ip = -^ + -(1 - p')ya{yfifl (37) 

This shows that r] plays the role of a volatility risk premium in that it enters as a drift adjustment for the 
volatility-driving process Y . However the nonlinearity of the PDE is through a quadratic term in tpy, as in 
the case of the entropic risk measure. 

Moreover, introducing the transformation 

'PQ{t,y) = — -log/(f,3;), 

y(l -p^) 

leads to the linear PDE problem for /: 

ft + Uy-ip + ilp')a(y)^(y)Y\f-\^\){l-p')f = 0, (38) 

fiT,y) = 1. (39) 

Therefore the indifference price is given by 

P{t, S,y) = ifit, S,y) + —r^-j- log f{t, y). 
7(1 -p^) 

In the numerical solutions, we take the volatility-driving process (Yt) to be an Ornstein-Uhlenbeck pro- 
cess with the dynamics: 

dYt - aim - Y,) dt + v ^Jl^ (p dWJ + p dWf), 

and we choose a function o-{Y) that gives realistic volatility characteristics. For the OU process (Yt), the rate 
of mean-reversion is a, the long-run mean-level is m and the long-run variance is v^. For the computations, 
we will take a = 5, m = 0, v^ = I, p = -0.2 and 

0.7 
o-{y) = — (arctan(y - 1) -i- n/2) + 0.03, 
t: 

so that o-{m) = 0.2050. The parameter v measures approximately the standard deviation of volatility fluctu- 
ations. The values are chosen such that the one standard deviation confidence interval for 7 is (-1, 1) and 
this translates roughly to the confidence interval (0.13,0.38) for volatility cr. The two standard deviation 
interval for volatihty is approximately (0.10,0.56). 
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We first solve the quasilinear PDE (ITTl) for ip using implicit finite-differences on the linear part, and 
explicit on the nonlinear part. Then we solve the linear PDE problem (l38l) for /. This procedure is done 
for a current stock price So = 100 and cr(Fo) = 0.223, calculating the solutions for European put options 
expiring in three months for strikes K in the range [70, 110] for different values of the distortion parameter 
77 and the risk-aversion parameter y > 0. 

Figure [3] reveals a more complex picture regarding the effect of r] away from the short maturity asymp- 
totic approximation. We see, as in Figure [T] from the asymptotics, increasing 77 increases the skew slope; 
however it also shifts down the levels of implied volatility around the money (as opposed to the opposite 
effect we saw in Figure [T]). 

Figure|4]shows, as we would expect, that increasing risk aversion y decreases the implied volatility skew 
which comes from the indifference price of the buyer who is willing to pay less for the risk of the option 
position. It also has a secondary effect of flattening the skew out of the money. 




0.1 -0-2 
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- - 11=0.3 

-^— T]=0.5 




Log Moneyness 



-0.4 -0.35 -0.3 -0.25 -0.2 -0.15 -0.1 -0.05 0.05 0.1 

Log Moneyness 



Figure 3: Implied volatility from the arctangent stochastic volatility model in terms of log-moneyness 
log(^/So) and 77 with fixed y - 0.5. 



4 Conclusion 



We have derived a nonlinear PDE for the implied volatility from indifference pricing with respect to dynamic 
convex risk measures defined by BSDEs under diffusion stochastic volatility models. Our asymptotic anal- 
ysis has highlighted the principal effect of the risk measure on option implied volatility at short maturities, 
namely through the appearance of f., in the first order correction solving (l27l) . 

In the example of Section 13.21 this translates explicitly to a steepening effect on the implied volatility 
smile from the distortion parameter 77. Numerical computations confirm this away from short maturity too, 
as well as quantifying the effect of risk aversion on the level of implied volatilities. 

The analysis can be used to infer some information about the driver, for example 77 and y in the family 
(|28] ). from market implied volatilities. This is illustrated using the short time asymptotics for the Hull-White 
model on S&P 500 impUed volatilities in Section [331 
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Figure 4: Implied volatility from the aixtangent stochastic volatility model in terms of log-moneyness 
log(.^/5o) and y with fixed 77 = 0.2. 



17 



5 Proofs 

Price- and Volatility Bounds 

Lemma 5.1 (First Price Comparison). Suppose that P = u- u and P - u - u for some u, u, Ti, u € 

C^'HQro,r) n C(e^), < TO < r, satisfy 

Pr>LP- Ug{-A(y),p'Ka(y)-Uy) - g{-A{y),p'Ka(y)Uy)), in Qr,,r (40) 



K 
1 

'k 



P,<LP- Ug{-A(y),p'Ka(y)u^ - g{-A{y),p'Ka{y)u;i), in Qr,,r (41) 



as well as 



Then P > Pon Q 



P>P on {0} X R^ n Qr^^r and ]0, tq] x R^ n <92ro,r- (42) 



ro,r- 



Proof. Even if the form of this comparison principle for sub- and superprices seems to be quite unusual, the 
proof follows directly along the lines of Friedman IIFri64[ Theorem 2.16] since the functions g contains no 
second derivatives. To be precise: this argument leads a version where the inequalities in (l40l ). (1421) and the 
conclusion are strict. But setting P - P + e(l + t) one gets a strict superprice and sending e to zero yields 
the stated version. n 

Proof of Proposition 12. lot 

Proof. To prove dTSl ). we intend to invoke the above comparison principle for the price process given in 
Theorem 12.91 since it is clear that the Black-Scholes pricing functions aie sub- resp. supersolutions of the 
PDE. Unfortunately we have the indifference price only as solution of a Dirichlet problem which does not 
give rise to directly comparable lateral boundary conditions, thus we have to alter the argument a bit. 
Denote for A'^ e N by u^-- the solution of the initial/boundary-value problems 



-u,'- + Lu"'^ = ^gi-A(y),p'Ka(y)u;^-y, 
u^'^0,x,y)= -{l-e-'y- 

""dBiO.N) 



u''-HT,x,y)\ = -P'''{T,x-a). 



a^' 



By a classical ai^gument IILSU67[ Section V.§8], we can extract a subsequence u^'"- of u^'- such that u^'"- 
converges together with its derivatives to u and it's derivatives pointwise in Lj. The same is true for U^ 
given by 

-Mf + LM^= ^gi-A(y),p'Ka(y)a^); 
a^iO,x,y)= 0; 

Thus P^'-(t, X, y) = u^'-{t, x, y) - u(t, x, y) satisfies 

-Pr'- + LP^-^ = ^gi-A(y),p'Ka(y)a^) - {g{-A(y),p'Ka(y)uf-); 

P^^^{Q,x,y)= (l-e'^)+; (43) 

^'''-(^'^'3')Lm^)= P''(j,x;a). 

and P^'- — > P along a subsequence. 

Noting that P^^ (t, x; a) is a subprice on every cylinder Qt,n by writing it in the odd form P^^ (t, x; a) = 

- (-P^^ (t, x; a)) to satisfy the comparison principle of Lemma [STTl we have P^^ (t, x;o;)< P^'- on Qt,n 



and hence in the Umit P^^ (r, x; a) < P. The other direction of Theorem 12.91 is proved, of course, by the 
same argument using P^^ (r, x; o^) as superprice. 

Finally a reformulation of the achieved result reads 

U{a;^T,x) < U{I{T,x,yfT,x) < U(a^T,x) 

and so the strict monotonicity of U with respect to the first variable yields ( |T9l ). n 

Deriving the PDE for the Implied Volatility 

To derive the PDE for the implied volatility, we note that from (IT3] ) and (IT4l ). it follows that the indiffer- 
ence price P obeys the equation 

-Pr + LP= j(g{-A(y),p'Ka(y)u,) - g{-A{y),p'Ka(y)uy)). (44) 

Substituting P{t, x,y) = U{I^{t, x,y)T, x) we note that the derivatives relate 

Pr=Ur- {tI\; P, = Ur ■ 2//,r + [/.,; Py - Ur ■ lIIyT- 

P^^ ^ Urr ■ AfllT^ + Ur ■ II^T + Ur ■ III^^T + Ur., ■ 4//,T + [/,,,,; 
P^_y ^ Urr ■ Al^hlyT^ + ^r ' llJyT + Ur ■ III^.T + Urx ' lUyT 
Pyy ^ Urr ■ AI^IyT^ + Ur ■ HyT + Ur ■ lUyyT. 

Plugging this into (l44l ) and dividing by Ur, we derive the equation (l24l) by noting that 

Uxx - Ux _ ^. UlL = ^ L_l. .^^_^ + l. _i_^2IIyi_ 

Ur ' Ur 2/4t2 2/t 8 ' Ur /t "^ 2 ' Ur My " "y ' 

The initial condition follows by observing that Vega vanishes as t tends to zero: 

Lemma 5.2. It holds that V - Uy — Uy — > uniformly on compacts as t ^ Q. 

Proof. Choose the cylinder Qj^r in the layer Lj such that the compact set is contained. Thus m, m € 
(-;'i+/?/2,2+/j^g^^-) Pi c{Lt) implies that Uy and Uy are yS/2-Holder continuous with some Holder constant c, 
whence 

\V{T,x,y)\ = \uyiT,y) - UyiT,x,y)\ 

< \uy(T,y) - Uy{0,y)\ + \uy{T,x,y) - Uy{0,x,y)\ < Ict^ -^ 

as T ^ since Uy{0, y) and Uy{0, x, y) exist and are equal to zero by the definition of the initial conditions, n 

Implied Volatility - Proof of the Main Theorem 

Lemma 5.3 (Second Price Comparison). Recall that u is the solution of the Cauchy problem (1141 ) and u of 
the same problem with initial condition equal to zero. Suppose that P, Pe C''2(Qro,r)nC(2To,r). < tq < T, 
satisfy 

Pr>LP- i[g(-A(y),p'Ka(y)Uy) - g{-A{y),p'Ka{y)Uy)), in Qr,,r (46) 



K 
1 

1 
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Lr^LP- Ug{-A(y),p'Ka{y)uy) - g{-A{y\p'Ka{y)uy)), in Qr,,r (47) 



as well as 

P>P>P on ]0, To] X r2 n dQro^r, (48) 

and 

P{0,x,y) - P{0,x,y) - P{0,x,y) = (1 - e^ 

then P> P> Pon Q^^.r- 



Proof. We note that inequality (1461 ) implies (P - P)r > L{P - P) which implies together with the lateral 
bound (P-P)>Oon ]0, tq] x R^ n dQrg^r and the initial condition P(0, x, j) - P(0, x, j) = that P > P on 
2ro,/- by the classical comparison principle. The second inequality is proved in the same way. n 



Lemma 5.4 (Volatility Comparison). Suppose that I, I_e C ' {Qtq,,-) n C(2ro,r). < tq < P, satisfy 
(tP), < t/M/ + P^/ + Tq{y)n - lTn\paiy)A(y) + -— 

^ • \ K{Uy - Uy) 



resp. 



/ 7^ ^ iKA-j 72^7 r <rj nTfl /^,/^ gi'^iyl p' Ka(y)uy) - g{-A(y), p' Ka(y)uy) 
(tI )r > tIMI + 1 Q1 + Tq{y)Uy - lTUy\pa{y)A{y) + 



K{Uy - Uy) 

in Qro,r together with the lateral comparison 

KT,x,y)<l{T,x,y)<l{T,x,y) on ]0, tq] X R^ n 5ero,r 
and the initial growth condition 



-2 



lim tI\t, X, y) = lim tI (t, x,y) = on {0} X R^ n 2^,,,^. (49) 

Then it holds that 

Ut, X, y) < 1{t, X, y) < 7(t, x, y) in Qr^,,- 

— -2 n 

Proof. Define first P(t, x,y) := U{I (t, x, 3')t, x) and P(t, x, j) := U(I_ (t, x, y)T, x). Then by the same calcu- 
lation as in the derivation of the PDE (E^ of the implied volatility in the paragraph above we get 

Pr>LP- j{gi-A(y),p'Ka(y)uy:) - gi-Aiy),p'Ka{y)uy)), in Qr,,r 

Lr^LP- j(g{-A(y),p'Ka{y)Uy) - g{-A(y),p'Ka(y)Uy)), in e,„,, 

as well as the lateral boundary condition 

P{t, X, y) < P{t, X, y) < P(t, x, y) in Qr,,r- 

Moreover, the growth condition ( |49l ) implies by the continuity of U that P(0, x, y) = P(0, x, y) = P(0, x, y) = 
(1 - e'^y. Thus we can use Lemma 1531 to infer P(t, x,j) < P{T,x,y) < P(T,x,y) in 2ro,r and the strict 
monotonicity of the function U in the first component yields the result. n 

Proof of Theorem l2l2t 
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Proof. Remark first that if tliere exists a solution to tlie PDE witfi some fixed initial condition, it has to be 
unique by the smoothness and strict monotonicity of U (and the boundedness of / proven in Proposition 
12.101 ) since otherwise the solution of the pricing PDE (Theorem |2.9b would not be unique. By the same 
reasoning we get also / € C^'^'^^^'^'^^iQr^r) ^ CiLj). Moreover, the eikonal equation (l25l) has a unique 
viscosity solution as proved in IIBBF04I Section 3.2]. To prove the theorem we will hence show that the 
solution of the eikonal equation is the only possible initial condition, i.e. that any solution of the PDE (|24l) 
has the eikonal equation as it's small time limit. More precisely we will show that there exist parametrized 
families of (time-independent) local super- and subsolutions of ((24l) which converge locally uniformly to 
the eikonal equation. This is done quite similar as in IIBBF04[ Section 3.4.], using an adapted vanishing 
viscosity method. However, in our setting the bounds on the volatiUty cr enable us to simplify the proof and 
circumvent some obscurities in the local volatility argument in IIBBF04L 

—E,6 

Define / ' (x, y) for e, 5 > as the solution of 

-6= -ir''f+i7'f^^+^A(in(r')y, 



\dB(m,r) ^^ 



r (50) 



where B{m, r) is an arbitrary disk. We will show that for r, 6, e > Q there exists a solution to this equation 
and for fixed r and 6 there exist eo > 0, tq > such that for all < e < eq this is a supersolution of (l24l ) in 

—E 6 

2ro,f Moreover, we show that I ' ^ f' locally uniformly as we send fist e and then 5 to zero. 
First step: Making the change of variables w := ln(/ ' ) one gets 

(f',^,^ = (r'f(.^o.)(p)>2p.,K.(^) j^)^ .„^o.)(p):) 

- o"^(3')(l - ^Wxf - 2po-(y)a{y)x{l - xwx)wy + a^(y)x^Wy =: Qw 

and the equation (l50l ) becomes 

-6 = -e^''' + §w + eAw; 
w\aB(m,r) = Ino". 

which admits a solution w e C^*^{B{m, r)) IILU681 Theorem 4.8.3] which is unique for sufficiently small r 
IILU681 Theorem 4.2.1]. 

— e S 

Second Step: By the Holder property of the derivatives of w resp. / ' implying the boundedness of the 
functions on the cylinder (as well as the boundedness of Uy and Uy and the differentiability of g) we can 
conclude that there exist constants C1-C4 solely depending on r such that 

1 O ,—E.&^ 7 -,— £.l5 

—y{I fgi < c,{r)T^ 



r/ Ml < C2{r)T 



r(q(y) 



g{-A(y),p'Ka{y)Uy) - g{-A{y),p'Ka{y)uy)\-s,s-s,s 



K{Uy - Uy) 



\-E,S-s,S , ^ 

Y ly <C3(r)T 
-£A(ln(7^''^))<£C4(r) 
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in B(m, r). We can conclude that 

—s,6 T. —£,6 n ,—£,(5,9 _ X . —s,6 . 

((T(/ )\={I f = 6 + 1(1 rg—+sAilniI )) 

>Tl MI +(/ ygz^-^Al fgi +Tq(y)I ly 

r, ,£5,£5^(-'^(3')'P'^«(3')"v) - gi-'iiyXp'Kaiy^y) 

— 2t1 ' 1' '■ '— 

K{Uy — Uy) 

+ 6 - (ci(r)T^ + C2{r)T + C3(r)T + £C4(r)j, 

—£,5 

thus for given 6 > and r > we can find indeed positive bounds on £o, tq sucli tliat / ' is a supersolution 
of (l24l ) for < e < eo in 2ro,r- In the same way one proves that we get for F'^ of 



S = -(ij^y + (f '^)^^^ + eA(ln (f '*)); 



V - l5B{m,r) — ■ 

subsolutions. 

Third Step: Having now super- and subsolutions, we can (by sake of Proposition 12. lOl ) invoke now the 
comparison principle Lemma [54] to conclude that 

fix, y) < I{T, X, y) < T\x, y) in Qr,,r 
for all < e < eo, eo and tq chosen as above. Thus we have 



J£,S/ 



-E,S. 



/ ix,y) < limini I{T,x,y) < lim sup /(r, ;c, j) < / ' {x,y). 

— £ 3 

Next we want to send £ to zero. Therefore we note that the families of solutions / ' , If'^ are uniformly 
bounded (by the constants V^ V o^ and V^ A o^ as a consequence of the comparison principle BGTOU 
Theorem 10.7. (i)] applied to the equations in the log-variables) and equicontinuous in e (since Holder- 

—E 6 

continuous with the same Holder constants). Thus by the Arzela-Ascoli theorem / ' converges along a 

subsequence uniformly on compacts to some limit function / e C°^^(2To,r)- This function is a viscosity 
solution of the PDE 

-6= -(ff^iffg^; 

\dB(m.r) ~ ^' 

An analogous result holds true for the subsolutions. Now sending 5^0, this gives by the same ai^gument a 
solution of the PDE 

which satisfies 7(0, y) = 0. Thus for t ^ 0, /(t, x, y) converges locally uniformly to I^ which is nothing else 
than the (by IIBBF041 Section 3.2] unique) viscosity solution of the eikonal equation (1251 ) with ij/ = x/I^. n 

References 

[ADEH99] P. Artzner, F. Delbaen, J. M. Eber, and D. Heath. Coherent measures of risk. Math. Finance, 
9:203-228, 1999. 

22 



[BBF04] H. Berestycki, J. Busca, and I. Florent. Computing the implied volatility in stochastic volatility 
models. Comm. Pure Appl Math., 57(10): 1352-1373, 2004. 

[BH08] P. Briand and Y. Hu. Quadratic BSDEs with convex generators and unbounded terminal condi- 
tions. Probab. Theory Related Fields, 141(3-4):543-567, 2008. 

[BK05] F. Benth and K. Karlsen. A PDF representation of the density of the minimal entropy martingale 
measure in stochastic volatility markets. Stochastics and Stochastics Reports, 77(2): 109-137, 
2005. 

[BK09] P. BaiTieu and N. Fl Karoui. Pricing, hedging, and designing derivatives with risk measures. 
In R. Carmona, editor. Indifference pricing. Theory and applications, pages 75-146. Princeton 
University Press, 2009. 

[CHOll] A. M. G. Cox, D. Hobson, and J. Obloj. Utility theory front to back - inferring utility from 
agents' choices. Preprint, 2011. 

[CHYP02] F. Coquet, Y. Hu, J. Memin Ying, and S. Peng. Filtration-consistent nonlinear expectations and 
related ^-expectations. Probab. Theory Related Fields, 123(1): 1-27, 2002. 

[DS06] F. Delbaen and W. Schachermayer. The mathematics of arbitrage. Springer Finance. Springer- 
Verlag, Berlin, 2006. 

[FFFIO] J. Feng, M. Forde, and J. -P. Fouque. Short maturity asymptotics for a fast mean-reverting 
Heston stochastic volatility model. SI AM J. Financial Math. , 1:126-141, 2010. 

[FJ09] M. Forde and A. Jacquier. Small-time asymptotics for implied volatility under the Heston 

model. Int. J. Theor Appl. Finance, 12(6):861-876, 2009. 

[FJLIO] M. Forde, A. Jacquier, and R. Lee. Small-time asymptotics for implied volatility under the 
Heston model: Part 2. Preprint, 2010. 

[Fri64] A. Friedman. Partial differential equations of parabolic type. Prentice-Hall Inc., Englewood 
Cliffs, N.J., 1964. 

[FS02] H. Follmer and A. Schied. Convex measures of risk and trading constraints. Finance Stock, 

6(4):429-447, 2002. 

[FSll] H. Follmer and A. Schied. Stochastic finance: An introduction in discrete time. Walter de 

Gruyter, 3rd edition, 2011. 

[Gia06] E. Rosazza Gianin. Risk measures via ^-expectations. Insurance: Mathematics and Economics, 
39(1): 19-34, 2006. 

[GTOl] D. Gilbarg and N. S. Trudinger. Elliptic partial differential equations of second order. Classics 
in Mathematics. Springer- Verlag, Berlin, 2001. Reprint of the 1998 edition. 

[HULOla] J.-B. Hiriart-Urruty and C. Lemarechal. Convex analysis and minimization algorithms II. 
Grundlehren der mathematischen Wissenschaften. Springer- Verlag, Berlin, 2001. 

[HULOlb] J.-B. Hiriart-UiTuty and C. Lemarechal. Fundamentals of convex analysis. Grundlehren Text 
Editions. Springer- Verlag, Berlin, 2001. Abridged version of Convex analysis and minimization 
algorithms I and II. 



23 



[KobOO] M. Kobylanski. Backward stochastic differential equations and partial differential equations 
with quadratic growth. Ann. Probah, 28(2):558-602, 2000. 

[KS05] S. Kloppel and M. Schweizer. Dynamic utility indifference valuation via convex risk measures. 
NCRR FINRISK working paper No.209. ETH Zurich, 2005. Shorter version published in Math. 
Finance 17(4):599-627, 2007. 

[LSU67] O. A. Ladyzhenskaya, V. A. Solonnikov, and N. N. Ural'tseva. Linear and quasilinear equations 
of parabolic type. Translated from the Russian by S. Smith. Translations of Mathematical 
Monographs, Vol. 23. American Mathematical Society, Providence, R.I., 1967. 

[LU68] O. A. Ladyzhenskaya and N. N. Ural'tseva. Linear and quasilinear elliptic equations. Trans- 
lated from the Russian by Scripta Technica, Inc. Translation editor: L. Ehrenpreis. Academic 
Press, New York, 1968. 

[MZ09] M. Musiela and T. Zariphopoulou. The spot and forward dynamic utilities and the associated 
pricing systems: Case study of the binomial model. In R. Carmona, editor. Indifference pricing. 
Theory and applications, pages 3^3. Princeton University Press, 2009. 

[Nut07] M. Nutz. Quadratic PDE and backward SDE. Diploma Thesis ETH Zurich, 2007. 

[Pen04] S. Peng. Nonlinear expectations, nonlinear evaluations and risk measures. Lecture Notes in 
Mathematics, 1856:165-253, 2004. 

[StalO] M. Stadje. Extending dynamic convex risk measures from discrete time to continuous time: A 
convergence approach. Insurance Math. Econom, 41 {3y.39l-404, 2010. 

[SZ05] R. Sircar and T. Zariphopoulou. Bounds & asymptotic approximations for utility prices when 
volatility is random. SI AM J. Control & Optimization, 43(4):1328-1353, 2005. 

[Tou07] A. Toussaint. Hedging and pricing with L^ convex risk measures in incomplete markets. PhD 
thesis, Princeton University, 2007. 



24 



